Temperature dependence of Coulomb drag between finite-length quantum wires 
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We evaluate the Coulomb drag current in two finite-length Tomonaga-Luttinger-liquid wires cou- 
pled by an electrostatic backscattering interaction. The drag current in one wire shows oscillations 
as a function of the bias voltage applied to the other wire, reflecting interferences of the plasmon 
standing waves in the interacting wires. In agreement with this picture, the amplitude of the current 
oscillations is reduced with increasing temperature. This is a clear signature of non-Fermi-liquid 
physics because for coupled Fermi liquids the drag resistance is always expected to increase as the 
temperature is raised. 
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Coulomb drag phenomena in coupled one- 
dimensional (ID) electron systems have been 
investigated quite extensively in the past 

0, S 1 H 0, S 0, 110. The interest has 
mainly been driven by the fact that Coulomb drag, 

1. e. the electrical response of one wire as a finite bias 
is applied to the other wire, seems to be an ideal 
testing ground for Tomonaga-Luttinger-liquid (TLL) 
physics in nature. This is because both inter-wire 
and intra-wire Coulomb interactions substantially 
modify transport properties such as the average 
current and the current noise. On the experimental 
side, there have been a few works, some of which 
have claimed to have observed TLL behavior in the 
drag data 11, 12, Recently, Yamamoto and 
coworkers have measured Coulomb drag in coupled 
quantum wires of different lengths and found peculiar 
transport properties that depend, for instance, on the 
asymmetry of the two wires LJ] . This experiment is 
the major motivation for our work. 

We analyze theoretically the Coulomb drag cur- 
rent of two electrostatically coupled quantum wires 
using the concept of the inhomogeneous TLL model 
HH, E3] ■ This model is known to capture the 
essential physics of an interacting ID wire of finite 
length coupled to non-interacting (Fermi liquid) elec- 
tron reservoirs. Within this framework, we are able 
to study finite-length and finite-temperature effects 
and therefore to make qualitative contact with the 
experimental setup of Ref. 14|. Since the Coulomb 



interaction varies between the wire regions and the 
lead regions, charge excitations feel the interaction 
difference at the boundaries and are known to exhibit 



Andreev-type reflections [171] . We show that these re- 
flections play a crucial role in the Coulomb drag setup 
illustrated in Fig. [T] Furthermore, we show that the 
quantum interference phenomena associated with the 
Andreev-type reflections considerably modify the drag 
current. This is particularly interesting as far as the 
temperature dependence of the drag current is con- 
cerned. For Fermi-liquid systems, it is well known 
that the drag resistance should always increase as the 
temperature is raised [6(. In our setup instead, the 
drag current at a fixed drive bias can either increase or 



decrease as a function of temperature. It crucially de- 
pends on the interference pattern due to finite-length 
effects. This is a clear signature of non-Fermi-liquid 
physics which could be observed in the double-wire 
setup of Ref. 0. 

The system considered consists of two in- 
teracting parallel wires (j = <,>) of fi- 
nite length L < (shorter wire) and L > (longer 
wire) connected to non-interacting semi-infinite ID 
leads (Fig. [T|) and is described by the Hamilto- 
nian H = Ej=<,> [ H j+ H Y] + H c- The intra- 
wire interaction is modelled through a TLL descrip- 
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with the piecewise constant interaction parame- 
ter gj(x) — g 3 < 1 in the wire region \x\ < 
Lj/2 and gj(x) — 1 in the non- interacting lead re- 
gions \x\ > Lj/2. The Fermi velocity vpj, the in- 
teraction strength gj, and the wire length Lj set the 
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FIG. 1: (color online). The system under consideration. 
Each interacting wire of length Lj [gray area, interaction 
parameter <7j(|a;| < Lj/2) = gj < 1] is connected to a pair 
of non-interacting leads > Lj/2) = 1]. The region 

of backscattering inter-wire interaction (red dashed box) 
extends over the length of the shorter wire. A voltage V is 
applied between the leads connected to the drive wire (here 
the longer wire j = >) and the backscattering-induced 
current in the drag wire (here the shorter wire j = <) 
is investigated. 
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FIG. 2: Drag current as a function of drive voltage for 
identical wires at zero temperature (solid curves). The 
inter-wire interaction strength ranges from strongly inter- 
acting (g = 0.1) to non- interacting (g = 1) for the dif- 
ferent curves (with = e\% s L 2 a ig /h 2 u>L and Vl = 
fiwi/e). The dashed curves show the dominant contri- 
bution oc V 49 ' 2 [given in Eq. JTOj] for g = 0.1,0.25,0.5. 



This term includes the contribution of the density^ 
density interaction which is most relevant to Coulomb 
drag [l|, Q when the Fermi wave- vectors of both wires 
are similar in magnitude, i.e. fcp< rs fcp> [Hi- 

In the following, we choose to apply a voltage V> — 
V to the longer wire (/ij* = — Mr = e V/2, drive wire) 
and none V < — to the shorter wire (/z^ = /x^ = 0, 
drag wire). The average current in the wires may then 

2 

be written as 7< = Idr and J> = \V — idr- In our 
model, the two currents i< and i> always flow in the 
same direction, which is due to momentum conserva- 
tion. This is known as positive Coulomb drag. 

In order to get an expression for the drag cur- 
rent Idr, we follow the formalism used in [lj| in the 
case of a single wire with an impurity. We consider 
the situation of weak inter-wire coupling. To second 
order in Abs, we obtain 

Idr = l£> / dR / drj dr (r,R), (5) 
Jo Jo 



frequency ujLj — i>fj/ 9jLj of the collective plasmonic 
excitations hosted in each wire. A voltage eVj — 
A*l "* — Mr" 1 applied to the leads is described by 



H 



dx jjLj(x)d x <&j(x), 



(2) 



with the piecewise constant electro-chemical potential 



Hj(x) = { 



for x < —Lj/2, 
for \x\ < Lj/2, (3) 
for x > Lj/2 



with fir^ 



A- 



This model is expected to cap- 
ture the essential physics of a quantum wire coupled 
smoothly to electron reservoirs (with typical smooth- 
ing length L s ) as long as L J=<!> 3> L s ^> Xp, where 
Xp is the electron Fermi wavelength [3, El- Fi- 
nally, we include an inter-wire backscattering inter- 
action over the length L < of the shorter wire, 

,L</2 

H c = A B s / da; cos{V47r[$ < (x)-$ > (x)]}. (4) 

J-L < /2 

I 



with the normalization ijj^ = eX% s L'^a^ < /H 2 ojl<, 
where a< = ujl</u c is the ratio between the plas- 
mon frequency of the shorter wire and a cutoff fre- 
quency lo c , of the order of the wire bandwidth. The 
plasmon frequency defines a voltage Vl = fiwL</e 
and a temperature Tl = hwL</k^. It is conve- 
nient to introduce corresponding dimensionless volt- 
age u = V/Vl and temperature 9 = T/Tl- The inte- 
grand in ||SJ), 



idr(r, R) 



2a%< 



x exp 



47rC<(r, R; r) + 4ttC* ; 



r R qr 



(6) 



involves the parameters I = L > /L < , p = g>/g<, 
and q = v-p > /v-p < . The correlation function Cj = 
Cj + Cj F of each wire can be decomposed in a 
zero-temperature and a finite-temperature contribu- 
tion given by fl8| 



Cf s (r,R;r) 
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fcez 
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■ , afc | / a j+ i(r- S r-2k) x | 2fe+I 

J V "3 + i(-2fe) ' 



ay + i(r - si? - 2A; - 1) s 
[ a 2 + (r-si?-2fc-l) 2 ] 1 /2^ 



(7) 



C"(r,i?;r) 



47T ^ 



| 2 fc| . / sinch[7r6l(T - sr - 2fc)] \ |2fe+i| , / sinch[7rg(r - sR - 2k - 1)] 
7j n ^ sinch[7r(9(-2fc)] ) lj n \smch[Tr6(r - sR- 2k - 1)] 



(8) 



with 7j = (1 — gj)/(l + gj) and sinchx = (sinhx)/x. 



It is to be noted that the expression for the drag cur- 



rent idr does not depend on whether the drive wire is 
the longer wire or the shorter one due to the symme- 
try of our model. In the following, we present results 
obtained by numerical evaluation of the triple integral 
involved in Q and (J6j> and discuss several analytical 
approximations. 

First we set the temperature to zero and consider 
identical wires (I = p — q = 1, thus we drop the wire 
index j). The drag current shows non-monotonous 
behavior and oscillations with period ~ l-Khio^ /e as a 
function of the bias voltage (Fig. [2|) . It decays at large 
voltages for g < 1/2 whereas it increases for g > 1/2. 
Thus, we obtain qualitatively the same behavior as in 
a dual Coulomb drag setup where a drive current is 
applied and a drag voltage is measured • Similar os- 
cillations as a function of voltage have been predicted 
in the context of two coupled fractional quantum Hall 
line junctions [201 ] . 

An analytic approximation can be derived in the 
limit u = eV/hujL 3> 1, that is for high voltages 
or long wires. In Eq. (|7|), the terms proportional 
to 7'™' account for contributions from plasmon exci- 
tations reflected \m\ times inside the wire. When the 
wire length is much longer than other relevant length 
scales, the contribution without any reflection m = 
becomes dominant and yields the integrand 



JdrO, R) 



7T 3 / 2 

u>i T(2g) 



(5) 



29-1/2 



■h g -i/2{ur), (9) 



where T(z) denotes the Gamma function and J„(z) 
the Bessel function of order v. The resulting expres- 
sion for the drag current, which involves hypergeo- 
metric functions, underestimates the amplitude of the 
oscillations with respect to the exact numerical re- 
sult. However, the behavior at large it, governed by 
the dominant contribution 



Jdr 



r(0) 
Z dr 



2r 2 (2 5 ) 



(?) 



49-2 



(10) 



shows good agreement in the appropriate parameter 
regime (dashed curves in Fig. Since we do per- 
turbation theory in Abs> the relation -C (e 2 /h)V 
has to hold. In the large it regime, this means 
(A B sV^c) 2 K/w c )(eI//fe i ; c ) 4 9- 3 < 1. 

Now we consider the situation where the two wires 
have different lengths. The qualitative behavior of the 
drag current does not change for increasing length ra- 
tio I = Ly/L <} but the peak positions get shifted to 
lower voltages (Fig. [3]). Here, neglecting plasmon re- 
flections in the correlation function of the wires leads 
again to the expression ([9]), which is independent of I. 
This fact explains why the drag current does not 
change appreciably as a function of I and indicates 
that the peak shifts observed result from plasmon re- 
flections inside the wires. Our studies are the first to 
analyze the effect of an asymmetry in the length on 
Coulomb drag phenomena in coupled quantum wires 
which is of recent experimental relevance [hH ]. 
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FIG. 3: (color online). Drag current as a function of the 
bias voltage and of the length ratio of the wires (with 
I&J = e ^ssL<a< /h 2 UL< and Vl = fko L </e). 



We now discuss in detail the temperature depen- 
dence of the drag current. For clarity, we consider 
again symmetric wires (I = 1) [23j]. The oscilla- 
tions of the drag current as a function of bias voltage 
get washed out with increasing temperature (Fig. [4]) . 
This behavior is consistent with the picture which at- 
tributes the oscillations to interferences of the plas- 
mon excitations of the wires. Thus, for bias voltages 
such that the drag current is close to a maximum at 
zero temperature, one observes a decrease of the drag 
current with increasing temperature, whereas the op- 
posite behavior can be observed close to a minimum of 
the zero-temperature drag current for strong interac- 
tions. This behavior is in stark contrast to the linear 
temperature dependence predicted for Coulomb drag 
between ID Fermi-liquid conductors [4j] and therefore 
a clear signature of TLL physics. Note that our Fig. |4| 
bears significant resemblance to Fig. 9 of Ref. 121 ]. 

A good approximation of the drag current can be 
obtained for temperatures much larger than the tem- 
perature associated with the plasmon frequency, 9 = 
k^T/fiWL 3> 1, by neglecting contributions from plas- 
mon reflections in the wires. Then, we obtain the 
dominant contribution 



J* ^ Jdr M fam H20J 4r 2 (2 5 ) 



(11) 



Taking the limit of low bias voltage in this result, 



Jdr 



dr 4r 2 (2 5 ) 



u, 



(12) 



we recover the power-law dependence T 4s ~ 3 of the lin- 
ear conductance predicted by renormalization group 
analysis {Ej]. At large bias voltage u > ^ > 1, we 
recover the zero-temperature result (|10p. 

In the case g = 1/2, Eq. (jTTJ) as well as the con- 
tribution to next order in (9 _1 can be brought into a 





FIG. 4: Temperature dependence of the drag current for 
identical wires with interaction strength g = 0.25. The 
solid curves (labelled a - e) are evaluated for tempera- 
tures given by T/Tl = kBT/huL = 0, 0.5, 1, 1.5, 5, respec- 
tively (with l£> = e\% s L 2 a 4g /H 2 u L and V L = fiwi/e). 
The dashed curve shows the high-temperature limit (|11[) 
for T/Tl = 5. The inset shows a similar plot for weaker 
intra-wire interaction g = 0.6 where the oscillations are 
less pronounced. 



FIG. 5: Drag current at high temperature (T/Tl = 5) 
for identical wires with interaction strength g — 0.5 (solid 

curve), where l£} = e\ BS L 2 a 49 /h 2 LOL and Vl = Swt/e. 
The dotted curve shows the dominant contribution in the 
high-temperature limit [first term in l|13[l]. the dashed 
curve includes the first correction as well [both terms 
in Q3)]. 



compact analytic form 
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where ^'(z) = d 2 lnr(z)/dz 2 denotes the trigamma 
function. The first term is Eq. (| 1 1 1) evaluated 
at g = 1/2 (dotted curve in Fig. O, and the second 
one is the dominant correction, which takes values 
within [—Ifa/9,0] (included in the dashed curve in 
Fig. This illustrates that the first-order approx- 
imation already yields a nice qualitative description 
for the full numerical result, which makes us confident 
that Eq. (| 1 1 1) is a good high-temperature approxima- 
tion also for g ^ 1/2. 

In summary, we have analyzed two coupled quan- 
tum wires that exhibit both a finite intra-wire interac- 
tion and a finite inter-wire interaction. We have taken 
into account finite-length effects within the inhomoge- 
neous TLL model that is known to capture the essen- 
tial physics of quantum wires coupled to Fermi-liquid 
reservoirs. We have investigated how an asymmetry 
in the lengths of the wires changes the drag current 
and we have predicted a rich temperature dependence 
of the drag current that shows clear signatures of non- 
Fermi-liquid physics. 

We would like to thank F. Dolcini, H. Grabert, 
M. Kindermann, Y. Nazarov, S. Tarucha, and M. Ya- 
mamoto for interesting discussions. This work 
was supported by the Swiss NSF and the NCCR 
Nanoscience. 



[1] K. Flensberg, Phys. Rev. Lett. 81, 184 (1998). 
[2] Y. V. Nazarov and D. V. Averin, Phys. Rev. Lett. 81, 
653 (1998). 

[3] A. Komnik and R. Egger, Phys. Rev. Lett. 80, 2881 
(1998). 

[4] V. L. Gurevich, V. B. Pevzner, and E. W. Fenton, J. 

Phys.: Condens. Matter 10, 2551 (1998). 
[5] R. Klesse and A. Stern, Phys. Rev. B 62, 16912 

(2000) . 

[6] V. V. Ponomarenko and D. V. Averin, Phys. Rev. 

Lett. 85, 4928 (2000). 
[7] B. Trauzettel, R. Egger, and H. Grabert, Phys. Rev. 

Lett. 88, 116401 (2002). 
[8] M. Pustilnik, E. G. Mishchenko, L. I. Glazman, and 

A. V. Andreev, Phys. Rev. Lett. 97, 126805 (2003). 
[9] G. A. Fiete, K. Le Hur, and L. Balents, Phys. Rev. B 
73, 165104 (2006). 
[10] M. Pustilnik, E. G. Mishchenko, and O. A. Starykh, 

Phys. Rev. Lett. 97, 246803 (2006). 
[11] P. Debray et al., Physica E 6, 694 (2000). 
[12] P. Debray et al, J. Phys.: Condens. Matter 13, 3389 

(2001) . 

[13] M. Yamamoto, M. Stopa, Y. Tokura, Y. Hirayama, 
and S. Tarucha, Physica E 12, 726 (2002). 

[14] M. Yamamoto, M. Stopa, Y. Tokura, Y. Hirayama, 
and S. Tarucha, Science 313, 204 (2006). 

[15] D. L. Maslov and M. Stone, Phys. Rev. B 52, R5539 
(1995). 

[16] V. V. Ponomarenko, Phys. Rev. B 52, R8666 (1995). 
[17] I. Safi and H. J. Schulz, Phys. Rev. B 52, R17040 
(1995). 

[18] F. Dolcini, B. Trauzettel, I. Safi, and H. Grabert, 
Phys. Rev. B 71, 165309 (2005). 

[19] The effect of a finite smoothing length L s on the con- 
ductance in an equivalent model has been analyzed by 
K. Janzen, V. Meden, and K. Schonhammer, Phys. 
Rev. B 74, 085301 (2006). 

[20] U. Ziilicke and E. Shimshoni, Phys. Rev. B 69, 085307 



(2004). 

A. Komnik and R. Egger, Eur. Phys. J. B 19, 271 
(2001). 

Inter-wire forward-scattering is not explicitely in- 
cluded in our model as it can be recast in a mere renor- 
malization of the intra- wire interaction strength g-j of 
each wire [jj. We also neglect electron tunneling be- 
tween the two wires for two reasons: (i) It is a less 
relevant process than He in a renormalization group 



sense [2lJ|. (ii) It can be tuned to zero in a Coulomb* 
drag experiment We assume to be away from 

half-filling where Umklapp scattering is forbidden. 
[23] In view of the results shown in Fig.|3] we do not expect 
any qualitative changes of the predicted temperature 
dependence for the case / = 1 as we make the wires 
asymmetric in length (I ^ 1). 



